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Abstract 

Wc present studies of double longitudinal spin asymmetries in semi-inclusive deep inelastic 
scattering using a new dedicated Monte Carlo generator, which includes quark intrinsic 
transverse momentum within the generalized parton model based on the fully differential 
cross section for the process. Additionally we employ Bessel- weighting to the MC events 
to extract transverse momentum dependent parton distribution functions and also discuss 
possible uncertainties due to kinematic correlation effects. 



1 Fully differential SIDIS cross section 

The study of the 3-dimensional structure of protons and neutrons is one of the central issues 
in hadron physics, with many dedicated experiments, either running (COMPASS at CERN, 
CLAS and Hall-A at JLab, STAR and PHENIX at RHIC), approved (JLab 12 GeV up- 
grade, COMPASS-II) or being planned (ENC/EIC Colliders). The transverse momentum 
dependent (TMD) partonic distribution (PDF) and fragmentation functions (FF) play a 
crucial role in gathering and interpreting information towards a true 3-dimensional imaging 
of the nucleons. TMDs can be accessed in several experiments, but the main source of infor- 
mation is semi-inclusive deep inelastic scattering (SIDIS) of polarized leptons off polarized 
nucleon. For SIDIS, the theoretical formalism is described in a series of papers [1, 2] using 
tree level factorization [3] where the standard momentum convolution integral [4] relates 
the quark intrinsic transverse momentum to the transverse momentum of the produced 
hadron PhT in semi-inclusive processes. 

In this work we present a model independent extraction of the ratio of polarized, (71, and 
umpolarized, /i, TMDs using a Monte Carlo (MC) based on fully differential cross section, 
in which we re-construct the final hadron transverse momentum after MC integration over 
the intrinsic quark transverse momenta. In the MC generator we used the model described 
in Ref. [1] that was numerically further evolved in Ref. [2]. The Bessel- weighted asymme- 
try, providing access to the ratio of Fourier transforms of gi and /i, has been extracted. 
The uncertainty of the extracted TMDs was estimated using different input models for 
distribution and fragmentation functions. 

A fully differential Monte-Carlo generator has been developed to describe the the SIDIS 
process when a final state hadron is detected with the final state lepton. 



l{l) + N{P)^l(l') + h{Ph)+X, 



(1) 



where i is the lepton, N the proton target and h the observed hadron (four-momenta 
notations are given in parentheses). The virtual photon momentum is defined q = I — V 
and its virtuahty = —q^. 

The fully differential SIDIS cross section used in MC is given by [1]: 



da ^ ^ 



dxdydzdp^dk^ 



^ J {fq{X: k±)Dq^h{z,P±) + Wl - egq{x, k±) Dq^h{z , P±)} 



(2) 



where the summation runs over the quark flavors and e, K(x, y) and J(a;, Q^, ki_) are some 
kinematic factors defined by the elementary scattering process [1], x is the Bjorken variable, 
/c_i_ is the initial quark transverse momentm, is the transverse momentum of the final 
hadron with respect to scattered quark, and y and z are the fractional energies of the virtual 
photon and detected hadron. For our studies we used simple factorized Gaussians for the 
fi{x,k±) and gi{x,k±) distribution functions and Di fragmentation function, with widths 
given by fits from available world data. 

Figure 1 shows the two dimensional plot of k^/Q"^ versus x. One can sec a clear cor- 
relation between the transverse and longitudinal momenta, where the dashed black curves 
[2] defining the upper bounds. Restrictions at large x come from the energy and momen- 
tum conservation and at small x from the requirement that the parton should move in 
the forward direction with respect to the parent hadron {kz > 0). The red curves (solid 
for g2 = 1 GeV^ and dashed = 3 GeV^) are calculated for a non-zero proton mass, 
requiring the parton to move in the forward direction with respect to the parent hadron (in 
Ref . [2] the proton mass was neglected) . 

We note that the portion of the data above the black dashed curve decreases with 
decreasing x. Thus, the smaller the x, the smaller is the difference between red solid and 
red dashed curves implying a diminishing role of the hadron mass at higher Q^. 



2 Bessel-weighted extraction of the double spin asym- 
metry All 

Extraction of the double spin asymmetry All, defined as the ratio of difference and sum 
of electroproduction cross sections for antiparallel, cr+, and parallel, , configurations 
of lepton and nucleon spins, has been performed using the Bessel- weighting procedure 
described in Ref. [5]. Within that approach one can extract the Fourier transform of 
the double spin asymmetry, A'^^'''^^'''^\bT), using measured double spin asymmetries as a 
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function of the PhT [6], for fixed x, y, and z bins. 

where 6t is the Fourier conjugate of the PhT- The Fourier transforms of hehcity dependent 
cross sections, a^{bT), can be extracted by integration (analytic models) or summation (for 
data and MC) over the hadronic transverse momentum, weighted by a Bessel function Jq, 

a^{bT)^S^ = J2jo{bTPhT,i). (4) 
i=l 

In the Fig. 2 the red points represent the outcome of the Bessel-weighted asymmetries 
from the MC sample, while the blue curve represents the analytical expression using 

< >g^ and <k']_>f^ from the fits to distributions from the same MC sample. 

Within the hx range of 6y < 5 ~ 6 GeV^^ ^ 1 fm the Bessel-weighted asymmetries 
could be extracted with a minimum of 2.5% accuracy, although with some systematic 
shift. The shift observed in the reconstructed value is due to the kinematical restriction 
introduced by energy and momentum conservation [2] , which deforms the Gaussian shapes 
of the k±_ and p±_ distributions. In experiment there is always a cutoff at high PhT due to 
the acceptance and low cross section, as well as a cutoff at small Phx^ where the azimuthal 
angles are not well defined. A correction factor, accounting for the missing P^t range 
above the maximum value accessible {Pmax) in a given experiment, was estimated, based 
on an analytic calculation of the contribution above that value. The systematic uncertainty 
of that correction could be estimated from the variation of the maximum Pmax within the 
resolution of the experiment. In Fig. 2 the blue filled squares represent this correction using 
the Gaussian distributions and the open black squares represent the numerical integration 
up to that exact A;^^^ that we used for correction. 
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Figure 1: The k\/Q'^ versus x = xb is presented for 11 GeV electron beam. Black dashed lines are from [2], with proton 
mass zero approximation. Events above black dashed curve are due to the hadron mass. Red curves are calculated under the 
fc^ > condition with non-zero hadron mass and for two values. Electron beam energy is 11 GeV. 
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Fi gUre 2: Extracted Bessel-weighted asymmetry versus bx with and w/o correction, compared to values calculated 
analytically and numerically directly from the input. The MC sample was produced assuming simple Gaussian DFs and EEs 
with < >9i= 0.8 < k']_>f-^ and gi{x) = /i(x)a;°'''. 
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